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We propose that FU Orionis outbursts may play an important role in 
maintaining the slow rotation of classical disk-accreting T Tauri stars. Current 
estimates for the frequency and duration of FU Orionis outbursts and the mass 
accretion rates of T Tauri and FU Orionis stars suggest that more mass may be 
accreted during the outbursts than during the T Tauri phases. If this is the case, 
then the outbursts should also dominate the accretion of angular momentum. 
During the outbursts, the accretion rate is so high that the magnetic field of 
the star should not disrupt the disk, and the disk will extend all the way in to 
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the stellar surface. Standard thin disk models then predict that the star should 
accrete large amounts of angular momentum, which will produce a secular 
^ <— | spinup of the star. 

We present boundary layer solutions for FU Orionis parameters which show 
that the rotation rate of the accreting material reaches a maximum value which 

H 

is far less than the Keplerian rotation rate at the stellar surface. This is due to 
the importance of radial pressure support at these high mass accretion rates. 
As a result, the rate at which the star accretes angular momentum from the 
disk drops rapidly as the stellar rotation rate increases. In fact, the angular 
momentum accretion rate drops below zero for stellar rotation rates which are 
always substantially below breakup, but depend on the mass accretion rate and 
on the adopted definition of the stellar radius. When the angular momentum 
accretion rate drops close to zero, the star will stop spinning up. Faster stellar 
rotation rates will produce negative angular momentum accretion which will 
spin the star down. Therefore, FU Orionis outbursts can keep the stellar 
rotation rate close to some equilibrium value for which the angular momentum 
accretion rate is small. We show that this equilibrium rotation rate may be 
similar to the observed rotation rates of T Tauri stars; thus we propose that FU 
Orionis outbursts may be responsible for the observed slow rotation of T Tauri 
stars. This mechanism is independent of whether the disk is disrupted by the 
stellar magnetic field during the T Tauri phase. 
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1. Introduction 

The FU Orionis stars are a class of accreting pre-main sequence stars which are 
observed to be T Tauri stars undergoing large outbursts (see Hartmann, Kenyon, & 
Hartigan 1993 for a review). During the outbursts, these systems brighten by ~ 5 — 6 
magnitudes over a period of ~ 1 — 10 years. They reach luminosities of a few hundred 
L , with spectra indicating maximum temperatures ~ 6500 — 7000 K. The outbursts last 
for ~ 100 years, although this number is poorly known, since no outbursting systems have 
been observed to return to their pre-outburst state. Hartmann & Kenyon (1985, 1987) 
showed that FU Orionis stars contain accretion disks by demonstrating that spectra of these 
systems contain double-peaked absorption lines. Simple steady accretion disk models of two 
FU Orionis systems, FU Orionis and V1057 Cygni, were constructed by Kenyon, Hartmann, 
& Hewett (1988, hereafter KHH). The models which best fitted the observed broad-band 
spectral energy distributions of these systems had accretion rates M ~ 1O~ 4 M yr _1 , and 
central stars with masses M* ~ 0.3 — 1 M and radii i?* ~ 4R . The presence of a disk 
in these systems suggests that their outbursts may result from a disk instability similar to 
those believed to cause similar outbursts in dwarf novae (Clarke, Lin, & Papaloizou 1989; 
Clarke, Lin, & Pringle 1990; Bell & Lin 1994; Bell et al. 1995). 

Disk accretion also has important implications for the rotational evolution of the 
accreting star. In the standard picture of disk accretion (Shakura & Sunyaev 1973), the 
star accretes angular momentum from the disk. The angular momentum accretion rate is 
simply the Keplerian specific angular momentum at the stellar surface multiplied by the 
mass accretion rate. This quickly spins up the accreting star. Thus, in the standard picture, 
if T Tauri or FU Orionis stars have accreted even a small fraction of their mass through a 
disk, they should be spinning rapidly. The rotation rates of FU Orionis stars are not well 
known, since the accretion luminosity is so high that it overwhelms the stellar luminosity. 
The rotation rates of T Tauri stars are more easily observable and have been studied in 
some detail (Bouvier et al. 1993, 1995). They tend to be about a tenth of breakup speed, 
on average, with rotation periods usually around 6-8 days. Since the angular momentum 
added by disk accretion should rapidly spin T Tauri stars up to breakup speed, some 
authors have hypothesized that the spinup of T Tauri stars is controlled by the interaction 
between the disk and a strong stellar magnetic field (Konigl 1991; Cameron & Campbell 
1993; Hartmann 1994; Shu et al. 1994). 
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This mechanism for maintaining the spin rate of a star near some equilibrium value 
was first explored in detail in order to explain the spin behavior of some accreting neutron 
stars (Ghosh, Lamb, & Pethick 1977; Ghosh & Lamb 1979a,b). In this model, the stellar 
magnetic field disrupts the disk at some distance from the stellar surface. The location 
of this disruption radius is determined by the strength of the magnetic field and the mass 
accretion rate; a stronger magnetic field disrupts the disk farther away from the star, while 
a larger mass accretion rate moves the disruption radius closer to the stellar surface. The 
Keplerian rotation speed at the disruption radius determines approximately how fast the 
star will rotate. In order to produce the rotation speeds of most T Tauri stars, the disk 
would have to be disrupted around 5 stellar radii from the stellar surface, and for the 
accretion rates commonly ascribed to T Tauri stars, this would require a stellar magnetic 
field ~ 1 kG (Kdnigl 1991). 

The magnetic hypothesis for explaining the slow rotation of T Tauri stars does not 
include the effects of angular momentum accretion during FU Orionis outbursts. Event 
statistics of these outbursts, in which the number of observed events is compared to the 
number of young stars in the observable volume, suggest that each star probably undergoes 
multiple outbursts (Herbig 1977; Hartmann & Kenyon 1985). Estimates for the time 
between outbursts vary widely from ~ 10 3 — 10 5 years (e.g. Bell & Lin 1994; Kenyon 1995). 
If the accretion rate during this period is the typical T Tauri rate of ~ 1CT 7 M yr™ 1 , then 
1CT 4 — 1CT 2 M Q will be accreted during the T Tauri phase. During the intervening FU 
Orionis outbursts, the accretion rate is ~ 1CT 4 M yr -1 for ~ 100 years, so ~ 10~ 2 M Q 
should be accreted during the outburst. This suggests that as an accreting pre- main 
sequence star evolves through multiple cycles of FU Orionis outbursts and T Tauri quiescent 
phases, most of the mass accreted by the star may be added during the outbursts. If this 
is the case, then we should expect that the outbursts will also dominate the accretion of 
angular momentum. 

During the FU Orionis outbursts, the stellar magnetic field should have little effect on 
the accretion disk. As noted above, in order to produce the observed rotation rates of T 
Tauri stars, the field should disrupt the disk at a radius of about 5 stellar radii. When such 
a star experiences an FU Orionis outburst, the mass accretion rate increases by a factor of 
order 1000. This should be sufficient to overwhelm the magnetic field, and allow the disk 
to reach all the way in to the stellar surface. In the standard picture of accretion disks, 
since the star accretes a large amount of mass during the FU Orionis outburst, it should 
also accrete a large amount of angular momentum, which will spin the star up substantially. 
This will in turn require that the star loses large amounts of angular momentum during the 
T Tauri phase in order to stay at the observed slow rotation rate. If the star accretes less 
mass in the T Tauri phase than it does in the FU Orionis outburst, it will be unable to lose 
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sufficient angular momentum to spin back down, and over time, the star will continue to 
spin up toward breakup speed. 

If the disk reaches the stellar surface, and the star is rotating slowly, a viscous boundary- 
layer forms between the rapidly rotating disk material and the star. The boundary layer 
controls the flow of angular momentum between the disk and the star. Also, if the central 
star is not rotating, the boundary layer should produce as much luminosity as the disk. 
Clearly, if we wish to understand the angular momentum transfer and to compare disk 
models to observed spectra of FU Orionis systems, it is essential to include the boundary 
layer in our disk models. 

Our initial investigation of the boundary layer structure examined the question of 
whether a star continues to accrete mass and angular momentum when it has been spun 
up to breakup speed (Popham & Narayan 1991, hereafter PN91). We found that as the 
star approaches the breakup rotation rate, the amount of angular momentum accreted 
per unit of mass accreted drops dramatically, and can even become negative, so that the 
star continues to accrete mass while losing angular momentum to the disk. Thus, there 
is some equilibrium rotation rate close to breakup for which the star accretes no angular 
momentum, so that it can continue to accrete mass without spinning up or down. This 
study used a very simple disk model with a polytropic relation between the disk pressure 
and density, so it was not clear that similar conclusions would hold when a more realistic 
treatment of the disk was used. More recently, we have explored the structure of the 
boundary layer region for several types of accreting stars, including T Tauri and FU Orionis 
stars (Popham et al. 1993, hereafter PNHK) and cataclysmic variables (Narayan & Popham 
1993; Popham & Narayan 1995), using a more realistic model which includes the energy 
balance and radiative transfer. 

In this paper, we examine the angular momentum transfer between the star and the 
disk for boundary layer solutions calculated for parameters corresponding to FU Orionis 
systems. The boundary layer model we use is similar to the one used in our studies of 
cataclysmic variables (Narayan & Popham 1993; Popham & Narayan 1995) and in our 
previous study of pre-main sequence stars (PNHK). We find that many of the conclusions 
reached by PN91 continue to apply when our more detailed model is applied to FU Orionis 
systems. When the stellar rotation rate is plotted against the height of the disk at the 
stellar surface, the solutions fall on multiple solution branches. We find that the angular 
momentum accretion rate drops as the stellar rotation rate increases, and that solutions 
exist for a wide range of angular momentum accretion rates, including negative values. 
There is an equilibrium stellar rotation rate for which the star can accrete mass without 
spinning up or down. The major difference from the PN91 results is that for these FU 
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Orionis-type solutions, the equilibrium stellar rotation rate is well below breakup. 

On the basis of these solutions, we propose that FU Orionis outbursts play an 
important role in the spin evolution of accreting pre-main sequence stars. The large amount 
of mass accreted during an outburst can change the stellar rotation rate substantially. If 
more mass is accreted during the FU Orionis outburst phases than during the intervening T 
Tauri phases, then the stellar rotation rate will stay close to the equilibrium rate established 
during the FU Orionis phase. Although this equilibrium rate is rather uncertain, we will 
show that for reasonable choices of parameters, it is comparable to the observed rotation 
rates of T Tauri stars. Thus, angular momentum transfer during the FU Orionis outbursts 
can explain the slow rotation of T Tauri stars without invoking strong stellar magnetic 
fields. If T Tauri stars have weak fields which allow the accretion disk to extend down to 
the surface of the star and spin up the star, negative angular momentum accretion during 
the FU Orionis outburst phases can keep the star spinning slowly. Alternatively, if T Tauri 
stars have strong magnetic fields which disrupt the disk, the FU Orionis phases will still 
dominate the spin evolution as long as they dominate the mass accretion. 

In §2, we discuss the model used to calculate our boundary layer and disk solutions 
for FU Orionis parameters, and the input parameters for our solutions. We describe a 
typical solution in detail in §3. We demonstrate the presence of multiple solution branches, 
and show how the character of the solutions varies along those branches. We also show a 
solution which accretes mass without accreting angular momentum. Finally, we show how 
variations in the mass accretion rate affect the solutions and solution branches. In §4, we 
describe the choice of an additional condition which allows us to find a relation between 
the stellar rotation rate and the angular momentum accretion rate. We discuss the spin 
evolution of T Tauri/FU Orionis systems and the implications of our results in §5. 



2. Disk and Boundary Layer Model 

2.1. Disk Equations 

The model we use to study the boundary layer uses a version of the "slim disk" 
equations developed by Paczyhski and collaborators (Paczyhski & Bisnovatyi-Kogan 1981; 
Muchotrzeb & Paczyhski 1982). This model is similar to the models used in previous papers 
(Narayan & Popham 1993; PNHK; Popham & Narayan 1995), so we will only describe it 
qualitatively here, and refer the reader to those papers for more details. 

The model assumes a steady state, so that the mass and angular momentum accretion 
rates M and J are constant with radius. We use the slim disk versions of the angular 



- 6- 



momentum, radial momentum, and energy equations. Thus, we dispense with some of the 
simplifying assumptions made in the standard thin disk treatment (Shakura & Sunyaev 
1973). We include pressure support and acceleration in the radial momentum equation, so 
that the rotational velocity Q is not required to equal the Keplerian value Qk- We allow 
the angular momentum accretion rate J to deviate from the standard value MQk{R,)R1 
by a factor j, so that J = jMf2^(i?Ji?^. Finally, we drop the assumption that the energy 
dissipated by viscosity is radiated locally at the same radius in the disk. Instead, we include 
terms for the radial transport of energy by radiation and by the accreting material itself. 
The radiative transfer scheme we use is the one described in PNHK. 

The modified angular momentum, radial momentum, and energy equations, along with 
the radiative transfer equations, provide a set of differential equations which can be solved 
with appropriate boundary conditions. The most important of these are the conditions on 
the angular velocity and the radial radiative flux. The angular velocity Q must match the 
stellar rotation rate fi* at R t , and the Keplerian value flxiRout) at R — R out = 100R t . At 
R t , there is an outward radiative flux characterized by an effective temperature T*, which 
we take to be 5000 K. This produces a luminosity of ~ 10 34 ergs s _1 entering the boundary 
layer from the star, which is a small fraction of the accretion luminosity L acc ~ 10 36 ergs s _1 . 
The equations are solved using a relaxation method. 

2.2. Solution Parameters 

Within the framework of the disk and boundary layer model described above and in 
the aforementioned papers, a large number of solutions can be found. These solutions are 
distinguished by various choices of input parameters. For our models, these include the 
mass accretion rate M, the angular momentum accretion rate J, the stellar mass, radius, 
and rotation rate M,, R t , and fi*, and the viscosity parameter a. In the case of FU Orionis 
systems, none of these parameters are directly observable, so they must be estimated by 
using the best current models to interpret the available observations. 

2.2.1. M,M t ,andR t 

These are the standard accretion disk parameters which determine the disk luminosity 
and temperature. Observed spectra of FU Orionis and V1057 Cygni were fitted with simple 
accretion disk models by KHH. Although these models did not include the boundary layer, 
they allow us to estimate the basic disk parameters as follows. First, the infrared flux is 
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produced primarily in regions of the disk which are fairly far from the stellar surface, so that 
the flux produced is insensitive to the value of R t . We have calculated simple blackbody 
spectra of our disk and boundary layer solutions for various choices of M and M s , and 
compared them to photometric data of V1057 Cygni. Fits to infrared photometry suggest 
that the combination MM t ~5x 1O~ 5 M 2 yr -1 . Our spectra and fits will be presented in 
a future paper (Popham et al. 1995). 

Another constraint is provided by spectra of FU Orionis and V1057 Cygni, which 
show that the maximum effective temperature reached in the inner disk is around 7000 K 
(KHH; Kenyon et al. 1989). In the standard Shakura & Sunyaev (1973) model, the disk 
reaches a maximum effective temperature T max = 0.287(GM t M/i?V) 1 / 4 at R = 49/36i?,. 
T m ax oc i?^ 3 ^ 4 ' so ft P rov ides a reasonable way of estimating R t ] however, the Shakura- 
Sunyaev formulation does not include any luminosity from the boundary layer, so the true 
temperature of the inner disk will be higher than T max . For MM t = 5 x 1O _5 M 2 yr" 1 , we 
find T m ax — 6500 K for R t = 3 x 10 11 cm = 4.31 R . This value for R t is somewhat larger 
than the observational estimates for T Tauri star radii, such as those tabulated by Bouvier 
et al. (1995), where R t ~ 1.5 - 2.5 R . 

This combination of parameters gives an accretion luminosity L acc ~ 
1.4 x 10 36 ergs s _1 ~ 35OL , which agrees with the bolometric luminosity estimates 
for most of the FU Orionis systems listed by Bell & Lin (1994). In the solutions presented 
in this paper, we keep M t = 0.5 M constant, and vary M and R t so that the accretion 
luminosity remains nearly constant. The three sets of solutions have M = 10~ 40 , 10 -4 ' 15 , 
and 1O" 4 - 3 M yr" 1 , and R f = 4 x 10 n ,3 x 10 11 , and 2.25 x 10 11 cm, respectively. 



2.2.2. a 

The viscosity parameter a is more difficult to determine. The most common way to 
do so is to try to estimate the viscous timescale in the disk by observing the disk outburst 
behavior. This assumes that the FU Orionis outbursts are produced by disk instabilities, 
or at least by some disturbance which dies away on the viscous timescale. Disk instability 
models of FU Orionis outbursts have found that small values of a ~ 10~ 3 — 10~ 4 were 
required to reproduce the long timescales of observed outbursts (Clarke, Lin & Pringle 
1990; Bell & Lin (1994); Bell et al. 1995). Unfortunately, these values of a lead to problems 
with the evolutionary timescales and gravitational stability of the disk (Bell et al. 1995). 
We have therefore adopted an intermediate value of a = 0.01 which we use for all of our 
models. 
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2.2.3. and j 

These parameters are usually ignored in accretion disk studies, since in the standard 
Shakura & Sunyaev-type model, the accreting star is assumed to be nonrotating, and the 
angular momentum accretion rate is just given by the Keplerian specific angular momentum 
at the stellar surface multiplied by the mass accretion rate. Our study treats the boundary 
layer region explicitly, so these assumptions are no longer necessary. We treat both the 
stellar rotation rate and the angular momentum accretion rate as free parameters for 
purposes of finding our solutions. Later, when we discuss the spin evolution of accreting 
pre-main-sequence stars in §4, we will attempt to define a relation which specifies the 
angular momentum accretion rate as a function of the stellar rotation rate. 

3. Solutions and Solution Branches 

3.1. A Typical Solution 

We begin by discussing an individual boundary layer solution with M = 10 -4 ' 15 M Q yr" 1 
and R t = 3 x 10 11 cm, with j = 0.90 and H(R t ) = 10 11 cm. These parameters are typical 
of most of the solutions presented here. Note that we have specified the disk height at R t 
rather than fi*; the reason for this is discussed below. Figure 1 shows the variation of the 
rotational velocity Q, the radial velocity vr, the effective temperature T e ff, and the disk 
height H in the region R = 1 — 3 R t . 

The rotational velocity f2 reaches a maximum around fl~5x 10 11 cm ~ 5/3 R t , so 
the dynamical boundary layer is about 2/3 R t wide. Note that the maximum Q is only 
about 70% of VL K at that radius. This means that even at the outer edge of the dynamical 
boundary layer, Vl 2 R is only about 50% of the gravitational term £l 2 K R; the other half of 
the support against gravity is provided by the pressure gradient. 

The radial velocity vr shows a similar profile to Q. It reaches a maximum of 
~ 3 x 10 4 cm s _1 at R ~ 6 x 10 11 cm ~ 2 R t . Thus, the accreting material moves radially 
through the boundary layer region in about 1 year. In the boundary layer, the radial 
velocity drops rapidly, reaching ~ —2000 cm s _1 at the stellar surface. 

The effective temperature T e // peaks closer to the stellar surface than either f2 or vr. 
In fact, the maximum T e // ~ 8300 K is reached at R = 3.8 x 10 11 cm, in the middle of the 
dynamical boundary layer. The radiative energy released in the dynamical boundary layer 
propagates out to larger radii, but most of it has been radiated from the disk surface inside 
2 R t . Thus, the thermal boundary layer (the region of elevated effective temperature) is 
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only slightly wider than the dynamical boundary layer. 

The disk height H is fairly large: the ratio H/R is close to 1/3 throughout the inner 
region. Note that in earlier studies of the boundary layer structure at lower accretion 
rates, we generally found that H increased rapidly in the pressure-gradient supported 
region inside the dynamical boundary layer. Here, as noted above, the accreting material 
is already substantially pressure-gradient supported at the outer edge of the dynamical 
boundary layer, and no rapid increase in H is seen at the innermost radii. 

3.2. Solution Branches 

One goal of this paper is to examine the variation in the boundary layer structure as 
the various system parameters discussed in §2.2 change. We are particularly interested in 
the spin history of accreting pre-main sequence stars, so we focus on the variations produced 
by changes in the rotation rate of the star and the angular momentum accretion rate j. 
Accordingly, we first varied while keeping all of the other parameters constant, including 
M, M t , R t , a, and j. This produced an intriguing result. We found that for some of our 
solutions, as we increased Q*, the disk height H t = H(RJ increased, while other solutions 
showed the opposite behavior: H t decreased as the star spun faster. 

We continued to systematically increase and decrease Q*, and plotted the values of H t 
in the resulting solutions as a function of fi*. This revealed that the tracks traced out in 
the fi* — H t plane connected together into a single curve, like the one shown in Figure 2. 
This curve has a Z shape; multiple solutions can be found for a single value of fi*. These 
solutions have different values of the disk height H t , and they are distinct solutions, despite 
the fact that they share the same values of all of the input parameters M, M t , R t , a, Q*, and 
j. In fact, by specifying H t rather than Q*, we can find a unique solution. The Z-shaped 
curve has three branches; on the middle branch, H t increases with increasing Q*, while on 
the upper and lower branches, H t decreases with increasing Q*. The existence of multiple 
solution branches for a single value of j agrees with the findings of our earlier work (PN91), 
which used a much simpler disk formulation. 

3.2.1. Variations Along the Branches 

The solutions along a Z-shaped curve show a wide range of variation in the structure 
of the boundary layer. In Figure 3 we have plotted 10 solutions with H m ranging from 
8 x 10 10 cm to 1.25 x 10 11 cm. The positions of these 10 solutions in the f2* — H m plane are 
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marked on the Z-shaped curve shown in Fig. 2 The solutions in the upper left have large 
values of H t and small values of fi* corresponding to a slowly rotating star. These solutions 
have a very broad boundary layer, as shown in Figure 3. The peak values of Q and vr sit 
out at R ~ 4i?„. The peak effective temperature, occurs much closer to the star, but the 
region of elevated effective temperature extends out to several stellar radii. The disk height 
ff is approximately proportional to R, with H/R ~ 0.4, as shown in Fig. 3c. 

Solutions with smaller values of ff t , but still lying on the upper branch, have 
progressively narrower boundary layers. The last solution shown on the upper branch, with 
if, = 1.05 x 10 11 cm, has VL max at R < 2R t (Fig. 3). Meanwhile, increases gradually, 
reaching a value of ~ 3.5 x 10~ 6 s _1 at the transition from the upper to the middle branch. 
Note that this value of fi* is less than 10% of Q K (R t ) ~ 5 x 10~ 5 s _1 . 

As H t continues to decrease along the middle branch, the boundary layer continues 
to shrink in radial extent. The changes in the solutions along the middle branch with 
decreasing if. seem to continue smoothly the changes along the upper branch, with one 
exception: fi* drops back down toward zero. The if . — 9 x 10 10 cm solution, near the 
transition from the middle to the lower branch, has Vt max at R ~ 4 x 10 11 cm ~ 4/3 R t (Fig. 
3), so the radial extent of the boundary layer has decreased substantially. This solution has 
O t ~5x 10~ 7 s- 1 ~ 0MQ K (RJ. 

Finally, along the lower branch, the trends noted above for the upper and middle 
branches continue. The boundary layer becomes quite narrow, with Q m ax at i? ~ 1.15ff„ in 
the H m — 8 x 10 10 cm solution. The value of fi* increases rapidly, and the solutions have a 
large angular velocity gradient dVt/dR at R t . This concentrated energy release produces a 
strong peak in T e ff near R t , but T e jj begins to drop as if. continues to decrease, because 
the large values of for these solutions means that a large portion of the rotational energy 
of the accreting material is retained rather than being dissipated in the boundary layer. 
The ratio H/R varies more with R for the lower branch solutions than for solutions on the 
other branches; the if. = 8 x 10 10 cm solution has H/R = 0.24 at R = R m but H/R ~ 1/3 
at R = 3R t . 

3.3. Angular Momentum Accretion Rate 

Thus far we have kept the angular momentum accretion rate j constant and showed 
that a family of solutions with a given value of j traces out a Z-shaped curve in the f2* — H f 
plane. Now we show the effects of varying the angular momentum accretion rate j. For each 
value of j, and using the same M, M t , R t , and a as before, we can follow the variations in 
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£7* as we vary H m . These variations are shown in Figure 4. 

As we increase j, we find a similar Z-shaped curve in the — H t plane, but the curve 
is displaced to the left, i.e. to smaller values of For j = 0.91, the transition between the 
lower and middle branches has moved to negative values of fi^. The middle-to- upper branch 
transition also moves left, but by a smaller amount. The Z-shaped curves for adjacent 
values of j nest inside one another, and they continue to be displaced to smaller fi^ as j 
continues to increase. 

If we decrease j, the curves are displaced to larger Q*, but they also begin to change 
shape dramatically. At j = 0.89, the middle branch has become much shorter and steeper, 
so that it extends over a very narrow range in fi*. At j = 0.84, the middle branch has 
vanished completely, leaving a smooth, monotonic curve with H t decreasing as increases. 

If we continue to decrease j, we continue to find solutions, even for j < 0. In these 
solutions, the star accretes mass while losing angular momentum. An example of a j = 
solution is shown in Figure 5. The angular velocity continues to rise all the way in to the 
stellar surface. Nonetheless, f2* is far below Qk{R*), so the inner region of the disk is 
almost completely pressure-supported. There is no maximum in Q, so in essence there is no 
boundary layer at all. The disk is fairly thick, with H/R ~ 0.4. We can find solutions for 
negative values of j, and even for large negative values such as j = —10. As j drops, the 
disk continues to thicken. We have included the curves for j = and j = — 1 in the fi* — H t 
plane in Figure 4. The curves are rather similar to the j = 0.84 curve discussed above; H t 
decreases monotonically with increasing VL*. 

3.4. Variations with M 

Having described the appearance of the family of curves in the Q* — H t plane produced 
by a range of values of j, we can examine the variations in this family of curves as the 
mass accretion rate changes. As discussed in §2.2, we also vary R t along with M so as to 
keep the accretion luminosity approximately constant, so that we examine 3 combinations 
of M and R,: M = 10~ 4 - , 10" 4 - 15 , and 1O" 4 - 3 M yr" 1 , and R t = 4 x 10 n ,3 x 10 11 , and 
2.25 x 10 11 cm, respectively. 

In Figure 6, we show the families of curves in the fi* — H t plane for M = 10~ 4 M Q yr -1 , 
i?„ = 4 x 10 11 cm and for M = 10~ 4 - 3 M yr -1 , R t = 2.25 x 10 11 cm. These can be compared 
to Fig. 4, which used intermediate values of M — 1O _4 ' 15 M yr -1 , R t = 3 x 10 11 cm. 
The shapes of the curves and their variation with j are quite similar for the 3 different 
accretion rates. The major difference is that the family of curves is displaced to larger 
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as M decreases, and to smaller fi* as M increases. Thus, at the lowest accretion rate, 
M = 1O~ 4 ' 3 M yr^ 1 , the transition from Z-shaped to monotonic curves occurs around 
0, ~ 1.5 x 1(T 5 s' 1 ~ 0.2Q K (RJ. This can be compared to the M = 10~ 4 - 15 M yr" 1 
case, where the same transition happens around fi* ~ 5 x 10 -6 s _1 ~ 0.1Qk{R*)- At 
M = 1O~ 4 '°M yr" 1 , all of the Z-shaped curves have moved to fi* < 0, corresponding to a 
star rotating in the opposite sense to the disk. 

4. Spin Evolution of Pre-Main-Sequence Stars 

4.1. Angular Momentum Accretion 

In order to study the spin evolution of a disk-accreting star, we need to know how much 
angular momentum is accreted by the star as a function of its rotation rate. The answer to 
this question will determine the spin evolution of the star. Of course, the spin evolution 
may be complicated by variations in the mass accretion rate or the mass and radius of the 
star. These considerations are important for T Tauri/FU Orionis stars, and we will discuss 
them later. 

To begin with, however, we concentrate on the variation of j as a function of f2* for a 
star with constant M, M t , and R t . The general sense of this variation is clear from Figs. 
4 and 6, which show how j varies with and H t . If we assume for the moment that H t 
stays constant, then j decreases as fi* increases. At small values of Q*, where the const ant-j 
curves are Z-shaped, j decreases gradually. As the star spins up, the constant- j curves 
become monotonic, and j decreases more rapidly. This suggests that the spin evolution of 
the star proceeds as follows: as long as j remains positive, the star continues to spin up. 
Eventually, j decreases to zero, and the star reaches an equilibrium rotation rate where it 
neither spins up nor spins down. 

The major uncertainty in this picture is in H t : for a given value of Q*, j also depends 
strongly on the value of H t . Thus, in some regions of the — H t plane, if H t were to 
decrease rapidly enough as fi* increases, j would increase. As we discussed in §3, we can 
find solutions by specifying j and either fi* or H t . In nature, one would expect that only 
the stellar rotation rate is an intrinsic parameter of the accreting system, and that 
both j and H t should be determined by the physics of the accretion flow. In order to 
reproduce this situation within our model, we need an additional condition which gives H t 
as a function of fi*. Such a condition will define a track in the f2* — H t plane, and the 
variation of j as a function of Q* along this track provides a basis for understanding the 
spin evolution of accreting stars. 
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4.2. Need for an Additional Boundary Condition 

The major reason for the uncertainty in H f is the lack of a clear definition of the stellar 
radius R t . In the solutions presented thus far, we have specified a value for R t , but we have 
not specified a boundary condition which ensures that the specified R t represents the true 
stellar radius in any physical sense. As a result, we can obtain a wide range of solutions 
which have very different physical characteristics at R„, as shown in Fig. 3. Thus, we 
need an additional boundary condition, a definition for R t which marks the location of the 
transition from the boundary layer to the star, based on the changes in the accretion flow. 

This is a difficult task, for several reasons. First, our model is set up in a way that 
specifically tries to avoid making distinctions between the disk, boundary layer and star. 
All of these regions are assumed to obey the same equations, with the same physical 
parameters. Second, the distinction between star and disk is much more difficult to make 
at large accretion rates. This can be seen clearly in Fig. 1. It is quite difficult to separate 
the accretion flow into different regions; Q, T e ff and H all vary smoothly over large length 
scales. This can be contrasted with the clear differences between the disk, boundary layer, 
and star which are visible in models for lower accretion rates, like those for cataclysmic 
variables (Popham & Narayan 1995) and T Tauri stars (PNHK). In those models, f2 
dropped rapidly in the boundary layer, and then stayed almost constant in the star; the 
transition between the two regions was quite rapid. Similarly, H began to rise rapidly with 
decreasing R inside the boundary layer. In the current solutions, the Q profile shows only 
weak indications of leveling off, and H continues to decrease approximately linearly with R. 
This makes it quite difficult to select a radius which represents a transition from boundary 
layer to star. 

Another perspective on this issue comes from examining the range of solutions in Fig. 
3, all of which use a single value of j. One way to locate the stellar radius is to select a 
solution from this set for which the boundary layer-star transition appears to occur at the 
inner edge of our computational grid, where R = R f . It seems clear that the solutions with 
the smallest values of H t are not appropriate. They have very large values of dVL/dR at R t , 
and it seems clear that if the solution were to continue in to R < R f , Vl would continue to 
drop precipitously. Thus, it appears that in these solutions, R f falls in the middle of the 
boundary layer, and the boundary layer-star transition would occur at R < R f . On the 
other hand, the solutions with the largest values of if„ have very broad regions, extending 
out to a few times R t , in which Q decreases quite gradually. It appears that these solutions 
have reached the boundary layer-star transition at R> R t , and the material then simply 
continues to settle as it flows in slowly to R f . The intermediate solutions with H t ~ 10 11 cm 
seem to have the best characteristics. Both fl and vr have dropped substantially but are 



-14- 



leveling off around R t . 

4.3. Choice of a Boundary Condition 

We have many possible choices for boundary conditions to specify R t . One possibility 
would be to try to match the shape of the disk surface to the shape of the surface of a star 
rotating at a speed however, as we have discussed, the disk height H does not vary in a 
way that resembles the transition from a disk to a star, so this would be difficult. The best 
variables to use for a boundary condition are probably the angular and radial velocities fl 
and vr. These variables show directly where the boundary layer is located, and compared 
to other variables like the central and effective temperatures of the disk, they depend far 
less on the details of radiative transfer, opacities, etc. 

Of the two velocities, we believe that the radial velocity provides a superior measure 
of the location of the boundary layer-star transition. This is because the angular velocity 
may or may not drop as the accreting material approaches the stellar surface; if the star is 
rapidly rotating, Q may drop only slightly, or may continue to increase, as we saw in the 
solutions with small and negative values of j in §3.3. In such cases, it becomes very difficult 
to distinguish the disk from the star on the basis of the Q profile. The radial velocity, on the 
other hand, always reaches a maximum value and then drops substantially as it approaches 
the stellar surface. Put another way, 0* does not have to be small, but vr^ does. 

The simplest condition we can put on vr is simply to select a value for vr ; *, and define 
R t as the radius where vr = vr^. This can then be used to determine the variation of 
if, and j with Q*. One simple way to do this is to keep j constant and vary H t , so as to 
move along one of the constant-j curves shown in Figs. 4 and 6. Fig. 3 shows that vr 
varies monotonically along a constant-j curve, such that vr decreases with increasing H t . 
Thus, by varying H f , we can find a solution for which vr = Vr^ at R = R f . Corresponding 
solutions can be found for other values of j, and the values of Vt* and H t for these solutions 
will then produce the desired relations for H t (Q*) and j(fi*). 

4.4. Variation of j with Vt* 

In order to illustrate how j varies with for this boundary condition, we choose a 
typical value of vr^ = 1000 cm s" 1 . This value is a factor of ~ 2 x 10 4 smaller than the 
free-fall velocity, and it is 5-15% of the maximum radial velocity reached in the resulting 
solutions. The variation of H for the solutions with vr = vr^ at R„ — 3 x 10 11 cm 
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is shown as a dashed line in Fig. 4. Note that H t is nearly constant for the full range 
of ft*, with H t ~ 1.07 — 1.11 x 10 11 cm, so that HjR t ~ 0.36. There is some variation 
of H m ; it decreases gradually from about 1.088 x 10 11 cm at ft* = to 1.075 x 10 11 cm 
at ft* ~ 1.3 x 10~ 5 s" 1 , and then rises more rapidly, reaching about 1.108 x 10 11 cm at 
ft* ~ 1.83 x 10~ 5 s -1 . 

More importantly, j decreases with increasing ft*; as the star spins up, it accretes less 
angular momentum. Figure 7a shows the variation of j as a function of ft*; j decreases 
gradually at first, and then drops more rapidly. The angular momentum accretion rate 
reaches zero at ft* ~ 1.5 x 10~ 5 s _1 , and rapidly drops to j = — 1 at ft* ~ 1.8 x 10 -5 s _1 . 
Note that these values of ft* are only about one third of the Keplerian angular velocity at 
this radius, £l K (Rj ~ 5 x 10~ 5 s _1 . Thus, the star stops accreting angular momentum 
when it is still rotating well below breakup. 

At other mass accretion rates, the variation of H t and j with increasing ft* is quite 
similar: H, stays approximately constant and j decreases gradually at first, then more 
rapidly. The major difference is that j decreases more rapidly at higher mass accretion 
rates, and reaches j = at a smaller value of ft*. This is illustrated in Figure 7b and 
c, where we plot for the M = 1O~ 4 -°M yr~ 4 and M = 1O~ 4 - 3 M yr" 1 solutions. 

We have used the same criterion, vr^ = —1000 cm s _1 , to define the stellar radius. At 
M = 1O _4 '°M yr _1 , j = for ft* ~ 6.6 x 10~ 6 s _1 , corresponding to a rotation period of 
~ 11 days. Also note that for this higher accretion rate, j is substantially smaller than one 
even at ft* = 0. At M — 10~ 4 - 3 M Q yr" 1 , ft*(j = 0) ~ 3 x 10~ 5 s" 1 , which corresponds 
to a much shorter rotation period of ~ 2.4 days. Recall that as M changes, we have also 
varied R t by a similar amount in order to keep the accretion luminosity constant. Thus 
the variation in ft*(j = 0) as a fraction of the breakup rotation rate ft^(i?J is not as 
dramatic as the variation in ft*(j = 0) itself; we find ft*(j = 0)/Q K (R m ) ~ 0.2,0.3,0.4 for 
M = 10" 4 - , 10" 4 - 15 , 1O~ 4 - 3 M yr" 1 , respectively. 

5. Discussion 

5.1. Spin Evolution of T Tauri/FU Orionis Stars 

We have found that boundary layer solutions for parameters corresponding to FU 
Orionis outbursts have small and negative angular momentum accretion rates for fairly low 
stellar rotation rates. This suggests that angular momentum transfer during the FU Orionis 
outbursts can keep the rotation rate of T Tauri/FU Orionis stars near to an equilibrium 
rotation rate ft* ie9 ,FC/ — ft*(j = 0). If the star is rotating faster than VL* m ,FU when it 
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experiences an FU Orionis outburst, it will lose angular momentum in the outburst phase 
and spin back down to £l* !eq ,FU- If, on the other hand, the star is rotating slowly when an 
outburst occurs, it will have j ~ 1 in the outburst phase and will spin up toward £l*, e q,FU- 
We have found that for accretion rates M = 5 x 10" 5 — 10~ 4 M Q yr -1 , j = for spin rates 
— 0) ~ 0.66 — 3.04 x 10~ 5 s -1 ~ Q*, e q,FU, corresponding to rotation periods of 2.4 — 11 
days. These are similar to the observed rotation periods of T Tauri stars (Bouvier et al. 
1993, 1995); however, = 0) is sensitive to the value assumed for v R ^, as discussed 
below. 

During the T Tauri phase, the stellar rotation rate will tend to move toward a different 
equilibrium value VL^ ^^t- If the stellar magnetic field is too weak to disrupt the disk, then 
the star will spin up with j ~ 1 until it reaches f2*, e(? ,TT — £Ik{R*)- If the stellar field does 
disrupt the disk, then Q*, e q,TT will be determined by the radius where the disruption occurs, 
and may be smaller or larger than VL*^ FU . Therefore j may be positive or negative, with 
a magnitude \j\ < (Rd/ R,) 1 ^ 2 (Ghosh, Lamb, & Pethick 1977), where Rd is the radius at 
which the field disrupts the disk. Models of magnetically disrupted disks in T Tauri stars 
suggest R d ~ 4 — 5R t , which would give j ~ ±2. 

The spin evolution of a T Tauri/FU Orionis star will depend on Q* t eq,FU, Q*, eq ,TT, and 
on the amounts of mass accreted in typical FU Orionis and T Tauri phases, AMpu and 
AMt T . If AMpjj > AMtt, as suggested by current estimates of the frequency and duration 
of FU Orionis outbursts, the stellar rotation rate will move toward f2* )eiJi TT during each 
T Tauri phase, but will return to £l* fiq ,FU during each FU Orionis phase. 

How large will the excursions in be? A rigidly rotating star has angular momentum 
J = /fi*, where / is the star's moment of inertia, so J = IQ* + iti*. If we assume that the 
moment of inertia changes slowly, so that J ~ /^X, then the star spins up and down at a 
rate Q*/QK{R t ) = (j /k)(M /M), where we have defined / = kMR 2 . Since k is substantially 
less than 1, the increase in as a fraction of the breakup rotation rate £Ik{R,) is several 
times larger than the fractional increase in the stellar mass M. If we assume j — 1, k — 0.2, 
then the accretion of AM = 0.01 M onto a 0.5 M star would increase fi* by 0.1Q K (R it ). 
If the typical change in during the T Tauri phase is larger than the difference between 
the equilibrium spin rates Q*, eq ,FU and fl*, eq ,TT, then f2* will reach f2* ie(? ,TT during each T 
Tauri phase, and remain there for the duration of the phase. 

It is important to note that even though the FU Orionis phases may produce more 
mass accretion than the T Tauri phases, AM FU > AM TT , a T Tauri/FU Orionis star 
will still spend most of its lifetime with moving away from Q*, e q,FU, due to the short 
duration of the FU Orionis phases. This means that the mean value of f2* during the T 
Tauri phase will be offset somewhat from fi* ieg ,F(7- The sign of this offset will depend on 
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whether the star spins up or down during the T Tauri phase, i.e. on whether Vt^ eqTT is 
smaller or larger than fl* t eq,FU- Also, the return to the T Tauri state after an FU Orionis 
outburst may be gradual, with the mass accretion rate declining slowly over a period of 
a few decades. At mass accretion rates intermediate between T Tauri and FU Orionis 
rates, the star should spin up, since the magnetic field may still be too weak to disrupt the 
disk, while the mechanism described in this paper will only operate at a large value of the 
equilibrium rotation rate. Despite these deviations, the high accretion rate of mass and 
angular momentum in the FU Orionis phase will return f2* to Sl*, e q,FU quite rapidly, and as 
long as AM FU > AM TT , £1* will continue to be close to Q* :eq ,FU- 

Note also that in general the equilibrium spin rate Jl* if >q will differ slightly from the 
value of £7* which gives j = 0. Since J = 7lX + we will have Cl* = for J = IQ*. 
This simply means that in equilibrium, J must be sufficient to compensate for changes in 
the star's moment of inertia. If the moment of inertia increases in response to accretion, 
j > at fi* >eg . If the moment of inertia decreases, as it might in an accreting white dwarf, 
then j(Q* teq ) < 0. For instance, if the star's radius of gyration stays constant, so that 
7/7 = M/M, then j(fl^ eq ) = kQ^ eq /Q K (Rj, where again 7 = kMR\ If Q^ eq /Q K (RJ is in 
the range 0.2-0.4 as suggested by our results, and k ~ 0.2, then j(£l* teq ) ~ 0.04 — 0.08, and 
Q*,eq will be only slightly smaller than = 0). 



5.2. Solutions with Negative Angular Momentum Accretion 

One of the most important results of this paper is that it establishes the existence of 
disk and boundary layer solutions with small and negative angular momentum accretion 
rates. In these solutions, the central star accretes mass but loses angular momentum. This 
is quite different from the standard thin disk formulation, in which the angular momentum 
accretion rate is always J = M£1 K (RJR^, or j = 1 in our units. The key to understanding 
this result is that the Q profiles of our solutions can be very different from the purely 
Keplerian fl assumed in the thin disk case. Most importantly, Q does not reach a maximum 
and then drop down to the stellar rotation rate Instead, Q continues to increase all the 
way in to the stellar surface at R = R t . The fact that Q is monotonic means that there is no 
radius Rn m ax at which dQ/dR = 0, where the viscous torque would vanish, and the angular 
momentum accretion rate would simply be the amount of angular momentum carried in 
the accreting material at that radius, J = MVL max R} lmax . Since for a steady flow J must be 
constant with radius, this would constrain J to be positive. In our low-j and negative-j 
solutions, there is no maximum in Q, so that the viscous torque carries angular momentum 
outward at all radii. If the torque exceeds the rate at which angular momentum is carried 
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in by the accreting material, the net angular momentum accretion rate will be negative. 

These solutions are very similar to the ones we found in an earlier paper (PN91). 
In that paper, we calculated the structure of the boundary layer using a simplified disk 
model with a polytropic pressure-density relation. We found that our solutions formed two 
branches in the fi* — H t plane for large values of j ~ 1. We did not find the lower branch, 
but seems likely that it exists for these solutions as well. The PN91 solution branches 
formed sharp corners where they met, unlike the rounded transitions between branches 
seen in the current solutions. This difference probably arises from the very different disk 
thicknesses in the two studies; the PN91 solutions generally had H/R < 0.01. We were also 
able to find solutions with small and negative values of j; for these values only a single 
branch was present in the f2* — if„ plane. Thus, despite the use of a different disk model, 
the qualitative behavior of these solutions is very close to that of the solutions described in 
this paper. 

One important question is why the solutions presented here reach negative angular 
momentum accretion rates when the stellar rotation rate is still well below the breakup 
rate ^(.RJ. For the PN91 solutions, we defined the stellar radius as the point where 
H/R = 0.1 in the accretion flow, and found that with this definition, j was very close to 1 
for most O*, but dropped precipitously when f2* ~ 0.915f2^(i?J. We have seen that the 
current solutions reach j = for ~ 0.2 — 0.5Qk{R„)- As discussed above, solutions with 
small or negative values of j must have dfl/dR < for all R. Thus, if Q ~ Qk close to R t , 
the star must be rotating close to breakup, f2* ~ Qk{R,), in order to allow a negative- j 
solution. If is much less than Qk(R«), ^ will have to decrease close to and there will 
be a maximum in Q, so that j — 1. 

In the FU Orionis solutions presented here, two factors combine to produce small and 
negative values of j when f2* <C f2^;(i?J. The first is the large radial width of the boundary 
layer. In the solutions with j ~ 1, Jl peaks fairly far from R t , so that the radial width of 
the boundary layer is comparable to R t . At the radius where Q peaks, Qk is substantially 
smaller than £l K (R t ); for instance, at R = 2R t , VL K ~ 0.35£l K (R t ). As the star spins up, 

only needs to reach this value to make dVL/dR < everywhere and allow j to be small or 
negative. Second, pressure support plays an important role, and Q is substantially smaller 
than Qk- This allows solutions with dfl/dR < everywhere for even smaller values of Q*. 
This also explains why we reach j = at smaller f2* for larger values of M; as M increases, 
pressure support plays a larger role, the disk becomes thicker, and fl is a smaller fraction of 
Vt K . 

It is also worth noting that the total accretion luminosity increases substantially as j 
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decreases. In general, the accretion luminosity varies according to the expression 



GMM 



J acc — 



R. 



(Popham & Narayan 1995). Thus, for large negative values of j, the accretion luminosity 
can be substantially larger than the standard value GM t M / R t . If the star spins up 
substantially during the T Tauri phase, then when an FU Orionis outbursts occurs, j 
may reach fairly large negative values. For instance, Fig. 4 shows that j = — 1 for 

~ 1.83 x 1CT 5 s _1 , which is only about a 6% increase in ^/^(.RJ over the point 
where j = 0. A 10% variation from Q^/Qk(R*) — 0.3 to 0.4, as discussed above, would give 
j ~ —2. This would give an accretion luminosity of ~ l.&&GM t M / R t) which would make 
the outburst even brighter than one would expect from the increase in M. This additional 
luminosity would come from the rotational energy of the star, which is released as the star 
spins down. 

Disk solutions with negative angular momentum accretion rates may be important in 
other types of accreting systems. In systems such as cataclysmic variables, the accretion 
rates are generally much smaller than those of FU Orionis systems. This produces a 
narrower boundary layer, and pressure support plays a smaller role, so we would expect 
that negative angular momentum accretion rates are only found for fi* ~ Qk(R„)- We are 
currently completing a study of negative- j solutions for disks around cataclysmic variables 
with M = 1O~ 8 M yr _1 which confirms these expectations (Popham 1995). In some other 
systems, such as embedded pre-main-sequence stars, or some X-ray binaries, the accretion 
rates may be large enough to produce solutions like the ones presented here, where negative 
angular momentum accretion rates are reached when the star is still spinning relatively 
slowly. 



5.3. Assumptions 

The results described above depend on a number of assumptions. One very important 
assumption is that FU Orionis outbursts occur in classical T Tauri stars; i.e., that the T 
Tauri and FU Orionis phases occur during the same epoch of pre-main sequence stellar 
evolution. Models of FU Orionis outbursts based on disk instabilities generally assume 
that the outbursts arise in T Tauri star disks, so that the FU Orionis and T Tauri phases 
correspond to the high and low states observed in other types of accreting stars, such as 
dwarf novae. The observational connection between T Tauri and FU Orionis stars is rather 
tenuous, since only one of the current FU Orionis systems was observed spectroscopically 
before it went into outburst, and no known FU Orionis systems have yet returned to their 
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pre-outburst levels. The one system for which a spectrum was taken prior to outburst, 
V1057 Cygni, resembled a typical T Tauri star (Herbig 1977). However, some of the 
observed characteristics of FU Orionis systems, such as their association with reflection 
nebulae and outflows, and their large far-infrared excesses, suggest that FU Orionis 
outbursts may arise in younger embedded sources rather than in classical T Tauri stars 
(Kenyon 1995). If this is the case, then the mechanism described here should serve to limit 
the stellar rotation rate during the embedded phase of pre-main-sequence stellar evolution. 
The rotation rate would then presumably be controlled by the stellar magnetic field during 
the T Tauri phase. 

We must also make some assumption about the location of the stellar radius, which is 
difficult to define in our model (see §4 for a discussion of this problem). We have chosen a 
particular value of vr^ to define R t , in order to provide an example of the variation of H t 
and j as a function of fi*. Unfortunately, the value of = 0) is fairly sensitive to the 
choice of vr : *. We have illustrated this in Fig. 4, where we have plotted lines corresponding 
to vr^ = —500 cm s" 1 and vr^ = —2000 cm s -1 . For these values of vr^, j reaches zero 
at f2* ~ 0.94 and 2.19 x 10~ 5 s _1 , respectively. This means that for a factor of 4 variation 
in vr^, = 0) varies by more than a factor of 2. Thus, even when we specify all the 
parameters of an FU Orionis system, we still have a substantial uncertainty in Q*(j = 0) 
resulting from the uncertainty in the definition of the stellar radius. As discussed in §4.5, 

— 0) also varies dramatically with M. Therefore we can only say that for what 
we believe to be reasonable choices of v R ^ and M, the resulting values of = 0) are 
comparable to T Tauri rotation rates. 

Our results are not strongly dependent on the magnetic field of the star, as long as the 
field is not strong enough to disrupt the accretion flow in the high-M FU Orionis phase. 
This would require a field substantially larger than the ~ 1 kG needed to disrupt the disk 
during the T Tauri phase. The disruption radius varies as Rd oc M" 2 ^ 7 B 4 ^ 7 (Pringle & 
Rees 1972), so to maintain the same Rd requires that B vary as M 1 / 2 . Since FU Orionis 
accretion rates are a factor of ~ 1000 higher than T Tauri rates, B would have to be about 
30 times larger to disrupt the disk, which would require 5~3xl0 4 G.A field this strong 
would disrupt T Tauri disks at R > 20i? t , which can probably be ruled out by observations. 

We have used stellar radii of 2,25, 3, and 4 xlO 11 cm, corresponding to 2.88, 4.31, 
and 5.75 R Q , in our calculations. These are larger than the radii generally found for T 
Tauri stars, which are in the range 1.5-2.5 R Q . These larger radii appear to be required by 
observations: smaller stellar radii would produce maximum temperatures in the inner disk 
which are substantially larger than those observed. This suggests that the accreting star 
may expand during the FU Orionis outburst. 
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Prialnik & Livio (1985) calculated the evolution of a fully convective 0.2 M star during 
the accretion of 2.5 x 1O~ 3 M of material at rates ranging from 10 _1 to 1O~ 1O M yr -1 . 
They also varied the fraction / of the accretion energy carried into the star with the 
accreting material; / ranged from 0.001 to 0.5. They found that for the accretion rates 
appropriate to FU Orionis outbursts, the star expanded stably if / < 0.03. If / > 0.1, 
the star expanded unstably. In our FU Orionis solutions, the central temperature of the 
disk at R t is ~ 1 — 3 x 10 5 K. The virial temperature for these parameters is ~ 10 6 K, 
so a reasonable fraction / ~ 0.1 — 0.3 of the accretion energy is advected into the star. 
This suggests that rapid stellar expansion may result from the high-M accretion during 
an FU Orionis outburst, but it is difficult to know how the mass of the star, the presence 
of rotation, and the departure from spherical symmetry would affect the Prialnik & Livio 
results. This expansion, and the subsequent contraction during the T Tauri phase, could 
affect the spin evolution of the star, and other aspects of the star's interaction with the 
disk, and clearly deserves further study. 

Our models also assume that the disk and boundary layer are in a steady state. This 
is clearly not true for FU Orionis systems, since they are experiencing outbursts. The 
viscous timescale is t V i SC ~ R 2 jv = a~ 1 (iJ/i?)~ 2 fi^ : 1 . For our solutions with a = 0.01, 
H/R ~ 0.3, this gives t visc ~ lOOOfi^ 1 . In the inner disk, we have Vl K ~ few x 10" 5 s _1 , so 
t V i SC ~ few x 10 7 s ~ 1 yr. This is comparable to the rise times of the fastest-rising systems, 
FU Ori and V1057 Cyg, and much shorter than their decline times. Thus FU Orionis 
systems should be reasonably well- modeled by steady state solutions for this value of a. 

5.4. Comparison with Observations 

Recent measurements of periodic photometric variations in T Tauri stars have clearly 
established that the classical disk-accreting T Tauri stars spin more slowly than weak-line 
T Tauri stars, which appear to lack accretion disks (Bouvier et al. 1993, 1995; Edwards et 
al. 1993). This has generally been interpreted to mean that the spinup of T Tauri stars 
is controlled by a stellar magnetic field strong enough to disrupt the disk. If FU Orionis 
outbursts regulate the spin evolution of T Tauri stars in the way described in this paper, 
it would also account for this observation. In both cases, disk accretion maintains the 
slow rotation speed of classical T Tauri stars. If FU Orionis outbursts cease late in the 
classical T Tauri stage, then the more rapid rotation of the weak-line T Tauri stars could 
be produced by accretion at the end of the classical T Tauri stage. 

Observations of the spin evolution of T Tauri stars could be used to test the ideas 
discussed in this paper. The most useful data would be measurements of the spin period 
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before and after an FU Orionis outburst. If the star spins down during an outburst, it would 
offer clear evidence that outbursts are regulating the rotation rates of T Tauri stars. If the 
star spins up, the result is less clear. If the amount of spinup is less than expected based on 
the amount of mass accreted during the outburst, it would suggest that the star is reaching 
the equilibrium rotation rate Q* t eq,FU at some point during the outburst. In this case the 
outbursts would still be limiting the rotation rate of the star. If, on the other hand, the 
star spins up as much as expected, it would suggest that the rotation rate must be limited 
by processes taking place during the T Tauri phase. One complication in interpreting the 
observations is the possibility that significant amounts of angular momentum might be 
carried off in the strong outflows which occur during outbursts. 

Observations of variations in the rotation rates of T Tauri stars which have not 
experienced observed outbursts could also be useful. Over a short timescale, variations 
in M could produce irregular period variations. Nonetheless, if FU Orionis outbursts are 
controlling the rotation rates of these stars, then their spin periods should vary secularly 
during the T Tauri phase as the rotation rate moves from £l* t eq,FU toward £l*, e q,TT- If 
instead the rotation rates of T Tauri stars remain constant over a long timescale, it would 
suggest that has stabilized at fi*, e <j,TT, which would indicate that the spin evolution is 
being controlled by angular momentum transfer in the T Tauri phase. 

The main sources of data which can be used to further constrain the solutions 
presented in this paper are the spectra of FU Orionis systems. These spectra yield valuable 
information on the disk temperatures and rotational velocities (Hartmann & Kenyon 1985, 
1987; KHH). The solutions presented in this paper were selected to have luminosities and 
temperatures in approximate agreement with those derived from observations. Because our 
main aim in this paper has been to understand the spinup and angular momentum transfer 
in FU Orionis systems, we have not attempted to match our solutions with observations 
in any detail. We defer this to a subsequent paper (Popham et al. 1995), where we plan 
to compare the spectra and rotational velocities produced by our disk and boundary layer 
solutions to those observed in FU Orionis systems. 
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Fig. I. — A boundary layer solution for typical FU Orionis parameters. The solution shown 
here has M = 10" 4 - 15 M yr"\ M, = 0.5 M , R m =3x 10 11 cm, j = 0.90, if, = 10 11 cm, and 
a = 10~ 2 . The four panels show the angular velocity Q, the radial velocity vr, the effective 
temperature T e ff, and the disk height if in the range from R = 3 — 9 x 10 11 cm = 1 — 3i? t . 
The dashed line shows the Keplerian angular velocity fix- 

Fig. 2. — Variation of the rotation rate f2* and the disk height if, at R — R t — 3 x 10 11 cm, 
for solutions with j = 0.90, M = 10~ 4 - 15 M yr" 1 , M, = 0.5 M , R m = 3 x 10 11 cm, and 
a = 10~ 2 . The solutions form a Z-shaped curve in the £7* — if„ plane, with three solution 
branches. 

Fig. 3. — Same as Fig. 1, but for the 10 solutions marked along the Z-shaped curve in Fig. 
2, with H t = 0.80,0.85,0.90, ...1.25 x 10 11 cm. The boundary layer width increases rapidly 
as if, increases, and angular velocity, radial velocity, and effective temperature all decrease. 

Fig. 4. — Same as Fig. 2, but for various values of j as labeled. All of the solutions on all 
curves have M = 10~ 4 - 15 M yr -1 , M, = 0.5 M , R t =3x 10 11 cm, and a = 10' 2 . As £7* and 
H t increase, j decreases, and the curves change from Z-shaped to monotonic. The dotted line 
marks the solutions with vr^ = —1000 cm s -1 , the upper dashed line vr^ = —500 cm s -1 , 
the lower dashed line vr,* = —2000 cm s -1 . 

Fig. 5. — A boundary layer solution with j — 0, where the star accretes mass but no angular 
momentum. This solution has M = 10~ 4 - 15 M yr" 1 , M = 0.5 M , R m = 3 x 10 11 cm, 
H t ~ 1.08 x 10 11 cm, and a = 10~ 2 . Note that v R ^ = -1000 cm s -1 , as discussed in §4.4; if 
this condition is used to define R t , then j = for ~ 1.55 x 10 -5 s -1 , corresponding to a 
rotation period of 4.7 days. 

Fig. 6. — Same as Fig. 4, but for two other choices of M and R t : (a) has M = 
1O" 4 -°M yr' 1 and R, = 4 x 10 11 cm = 5.75 R , and (b) has M = 10- 4 - 3 M o yr" 1 and 
R t = 2.25 x 10 11 cm = 3.23 R . Note that as M increases, the Z-shape of the curves 
becomes less pronounced and only appears at lower values of The dotted lines mark the 
solutions with vr^ = —1000 cm s _1 . 

Fig. 7. — Variation of the angular momentum accretion rate j with the stellar rotation rate 
£7*, for a set of solutions with a given value of v R at R t : (a) is for M = 10~ 4 - 15 M yr -1 , 
R t = 3 x 10 11 cm, and shows for v R ^ = -500, -1000, and -2000 cm s _1 ; (b) is 

for M = 1O~ 4 -°M yr" 1 , R t = 4 x 10 11 cm, and v R} * = -1000 cm s _1 ; and (c) is for 
M = 1O' 4 - 3 M yr -1 , R„ = 2.25 x 10 11 cm, and v R ^ = -1000 cm s" 1 . In all cases, j drops 
with increasing f2* and reaches j = for £7* well below £1k(R*)- 
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